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Abstract 

In this paper, we investigate the Rayleigh- Taylor instability problem for two compressible, 
immiscible, inviscid flows rotating with an constant angular velocity, and evolving with a free 
interface in the presence of a uniform gravitational field. First we construct the Rayleigh- Taylor 
steady-state solutions with a denser fluid lying above the free interface with the second fluid, then 
we turn to an analysis of the equations obtained from linearization around such a steady state. 
In the presence of uniform rotation, there is no natural variational framework for constructing 
growing mode solutions to the linearized problem. Using the general method of studying a family 
of modified variational problems introduced in [4] , we construct normal mode solutions that grow 

exponentially in time with rate like e*"^ c '^' _1 , where £ is the spatial frequency of the normal 
mode and the constant c depends on some physical parameters of the two layer fluids. A Fourier 
synthesis of these normal mode solutions allows us to construct solutions that grow arbitrarily 
quickly in the Sobolev space H k , and lead to an ill-posedness result for the linearized problem. 
Moreover, from the analysis we see that rotation diminishes the growth of instability. Using the 
pathological solutions, we then demonstrate the ill-posedness for the original non-linear problem 
in some sense. 

Keywords and Phrases: Rayleigh- Taylor instability; rotation; Hadamard sense. 
AMS Subject Classifications: 35L65, 35L60. 



1 Introduction 

Hydrodynamic instabilities at the interface of two materials of different densities are a critical issue 
in high energy density physics (HEDP). The Rayleigh- Taylor instability (RTI) occurs when a fluid 
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accelerates another fluid of high density [lOj [TTJ [13]. The RTI is ubiquitous in HEDP, such as 
high Mach number shocks and jets, radiative blast waves and radioactively driven molecular clouds, 
gamma-ray bursts and accreting black holes, etc. Due to the importance in physics mentioned above, 
there have been many studies related to RTI from both physical and numerical simulation points 
of view in the literature, but only few analytical results. One classical case of RTI is to consider 
two completely plane-parallel layers of immiscible inviscid fluids, the heavier on top of the light one 
and both subject to the earth's gravity. The equilibrium here is unstable to certain perturbations or 
disturbances. For this model, Y. Guo and I. Tice [3] established a variational framework for nonlinear 
instabilities. They constructed solutions that grow arbitrarily quickly in the Sobolev space by the 
method of Fourier synthesis, which leads to an ill-posedness to the perturbed problem. Ever since 
then, several works studied the influences of different physical quantities to the linear RTI, such as 
the stabilized effect of viscosity and surface tension (see [HE]), and the effect of magnetic field (see 
m E!), etc. 

While, if the two fluids are all subject to a uniform rotation with an constant angular velocity, 
how does this rotation influent the RTI? Stabilize it or make it more unstable? To the best of our 
knowledge, physicists [TJ showed that rotation can only slow down the growth rate of the disturbance 
a little bit, but not prevent the incompressible fluids from becoming unstable. Now, an interesting 
question is whether it is still true for compressible flows? In this paper, for inviscid compressible flows 
without the centrifugal force, we first prove that the linearized system is unstable in the Hadamard 
sense. Then, we establish the ill-posedness for the original non-linear problem in some sense. 

Next, we formulate the problem in details for further discussion. 

1.1 Formulation in Eulerian coordinates 

We suppose that the fluids are confined between two rigid planes. As in [3], we denote this infinite 
slab by Q, = M 2 x (— m, I) C M 3 . The fluids are separated by a moving free interface £(i) (t > 0) that 
extends to infinity in every horizontal direction. The interface divides f2 into two time-dependent, 
disjoint, open subsets £l±(t), so that f2 = f2+(t)UJ7_(t)U£(i) and £(i) = 0+(i)nf)_(t). The motion 
of the fluids is driven by the constant gravitational field along e3-the #3 direction, g = (0, 0, —g) with 
g > and the rotation with an angular velocity u = (0, 0,u>) about the vertical direction. Quantities 
describing fluids are their density and velocity, which are given for each t > by, respectively, 

p±(t, •) : n±(t) -> M + , u±(t, •) : Q±(t) -»■ K 3 . 

We shall assume that at a given time t > these functions have well-defined traces onto ~S(t). 
The fluids are governed by the following equations: 

d t p± + div{p±u±) = 0, 

p±{d t u± + u± • Vw±) + 2p±(u x u±) + Vp(p±) = -gp±e 3 , 

for t > and x £ Q±(t). Here we have written g > for the gravitational constant, e% = (0,0, 1) 
for the vertical unit vector, and — ge% for the acceleration due to gravity, 2p±(u x u±) represents the 
Coriolis force, while the centrifugal force p±X7\u> x x| 2 /2, like in (9j [15] , is neglected on the basis of a 
standard argument which indicates that the motion is dominated by gravitation and the centrifugal 
force is small. 

We assume a general barotropic pressure law of the form p± = p±(p) > with p± E C°°(0, 00) 
and strictly increasing. We will also assume that l/p'± 6 L^ c (0, 00). Finally, in order to construct a 
steady state solution with an upper fluid of greater density at S(t), we will assume that 

Z := {z £ (0,oo) I p-{z) > p+(z) and P-(z) G p+(0,oo)} / 0. 
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In particular, this requires that the pressure laws be distinct, i.e. p- ^ p + . 

Now, we prescribe the jump conditions that, from the physical point of view, both normal compo- 
nent of the velocity and the pressure are continuous across the free interface (since no surface tension 
is taken into account), see [HE). Therefore, the jump conditions at the free interface read as 

KHIs(t)=0, [P]|s(t)=0 

for each t > 0, where we have denoted the normal vector to by v, the trace of a quantity / on 
E(£) by /Ism and the interfacial jump by 

[/]|s(t) : = /+|s(t) - /-|s(t)- 

On the fixed boundaries, we consider that the normal component of the fluid velocity vanishes, 
that is, 

u + (t, x , — m) • e3 = U-(t, x , I) ■ e$ = 0, for all t > 0, x := (xi, X2) G M 2 . 

The motion of the free interface is coupled to the evolution equations for the fluids (jl.ip by 
requiring that the interface be advected with the fluids. This means that the velocity of the interfcae 
is given by (u • u)v. Since the normal component of the velocity is continuous across the surface there 
is no ambiguity in writing u ■ v. The tangential components of u± need not be continuous across 
and indeed there may be jumps in them. This allows possible slipping: the upper and lower 
fluids moving in different directions tangent to Since only the normal component of the velocity 
vanishes at the fixed upper and lower boundaries, {x 3 = 1} and {23 = — m}, the fluids may also slip 
along the fixed boundaries. 

To complete the statement of the problem, we must specify initial conditions. We give the initial 
interface S(0) = So, which yields the open sets ^±(0) on which we specify the initial data for the 
density and the velocity 

(P±,«±)(0,0 :ft±(0)^(R + ,M 3 ). 
To simply the equations we introduce the indicator function x an d denote 

P = P+Xn+ + P-Xn- , u = u + xn+ + , P = P+Xn+ + P-Xfi- ■ 

Hence the equations are replaced by 

dtp + div(pu) = 0, ^ 
p(d t u + u -Vu) + Vp(p) = -gpe 3 + paju 2 ei - puuie2, 

for t > and x £ f2/£(£). It will be convenient in our subsequent analysis to rewrite the momentum 
equations in (II. 2p by using the enthalpy function 

h(z) = f Z P -^dr. 
Ji r 

The properties of p ensure that h 6 C°°(0, 00). Thus, (|1.2p can be rewritten as 

d t p + divu = 0, 

dtu + u ■ Vu + V/i(p) = —g&z + 2uju2&\ — 2uju\e2- 
In the subsequent analysis it will be convenient to use the notation 

If} := f+\x s =0 - f-\x 3 =0- 
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1.2 Steady-state solution 

In order to produce RTI, we first seek a steady-state solution with u± = and the interface given 
by {x 3 = 0} for all t > 0. Then 0+ = := M 2 x (0,/), = =: M 2 x (-m,0) for all t > 0, 

and the equations reduce to the ODE 

= ~gp± m fi±, (1.4) 



subject to the jump condition 



dx 3 



P+(P+) = P-(P-) on {^3 = 0}. (1.5) 



Such a solution depends only on x$, so we consolidate notation by assuming that p± are the restric- 
tions to (0,/) and (— m, 0) of a single function po = po(xs) that is smooth on (— m,0) and (0,/) with 
a jump discontinuity across {X3 = 0}. 

From the assumption of the pressure function p and the definition of the set Z, there exist two 
positive constants I and m and a solution po to ()1.4|) f|1.5[) such that 

• po is bounded above and below by positive constants on (— m, I), and po is smooth when 
restricted to (— m, 0) or (0, /). 

• [PoJ = Po U=o - Po U=o > 0. 

Please refer to Section 1.2] for more details concerning construction of such a solution. In this 
paper, we always assume that I, m and the solution po satisfy the above properties. 

1.3 Formulation in Lagrangian coordinates 

Since the movement of the free interface ~S(t) and the subsequent change of the domains Q±(t) 
in Eulerian coordinates result in mathematical difficulties, we switch our analysis to Lagrangian 
coordinates as usual, so that the interface and the domains stay fixed in time. To this end, we define 
the fixed Lagrangian domains = M 2 x (0, 1) and f2_ = M 2 x (— m, 0). We assume that there exists 
invertible mappings 

77° :O±^O±(0), 

so that £0 = 7/^({x3 = 0}), {x3 = 1} = t/J_({x3 = /}), and {X3 = — m} = rp_ ({x^ = — m}). The first 
condition means that S is parameterized by either of the mappings rf^ restricted to {x 3 = 0}, and 
the latter two conditions mean that r{± map the fixed upper and lower boundaries into themselves. 
Define the flow maps r/± by the solution to 

d t r]±{t,x) = u±(t,r]±(t,x)), 
7]±(0,X) = (x). 

We think of the Eulerian coordinates as (t, y) with y = r](t,x), whereas we think of Lagrangian 
coordinates as the fixed (i, x) € M + x f2, this implies that £l±(t) = r/±(t,Q±) and that E(i) = 
r/ + (t, {2:3 = 0}), i.e., that the Eulerian domains of upper and lower fluids are the image of Q± under 
the mapping rj± and that the free interface is parameterized by r/ + (t,-) restricted to M. 2 x {0}. In 
order to switch back and forth from Lagrangian to Eulerian coordinates we assume that r]±(t,-) is 
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invertible. Since the upper and lower fluids may slip across one another, we must introduce the slip 
map S± : R+ x R 2 -> R 2 x {0} C R 2 x (-m, /) defined by 

S_(t,x') = r]Z\t,7] + (t,x',0)), x € IR 2 (1.6) 

and S+(t, ■) = SZ 1 (t,-). The slip map S- gives the particle in the lower fluid that is in contact with 
the particle of the upper fluid at x = (xi,x 2 ,0) on the contact surface at time t. 
Setting rj = X+V+ + X-V-j we define the Lagrangian unknowns 

(v, q){t, x) = (u, p)(t, V (t, x)), (t, x) € R + x (Q/{x 3 = 0}). 

Defining the matrix A := (Aij) 3x3 via A T = (Dr?) -1 := {djrii)^ 3 , and the identity matrix I = 
(Iij)sx3, then in Lagrangian coordinates the evolution equations for rj, v, q are, writing dj = d/d Xj , 

' d t f] = v, 

< d t q + qtr(ADv) = 0, (1.7) 
, d t v + AVh(q) = -gAVr]3 + 2uv 2 AVrn - 2wwiAVr/ 2 - 

Since the boundary jump conditions in Eulerian coordinates are phrased in terms of jumps across 
the interface, the slip map must be employed in Lagrangian coordinates. The jump conditions in 
Lagrangian coordinates are 

f (v + (t,x',0) -v-(S-(t,x',0))) -n(t,x',0) = 0, 
\p(q + {t,x',0)) =p(q-(S-(t,x',0))), 

where we have written n = v o r), i.e., 

dir)+ x d 2 7] + 



n := 



\dir] + x d 2 f]+\ 

for the normal to the surface £(i) = r] + (t, {x^ = 0}). Note that we could just as well have phrased 
the jump conditions in terms of the slip map S+ and defined the interface and its normal vector in 
terms of ry_ . Finally, we require 

V- (t, x' , —m) • es = v+(t, x f , I) ■ = 0. 

Note that since dtf] = v, 

e 3 ■ r] + (t,x',l) = e 3 ■ 77° (x',l) + / e 3 • v+(t, x, l)ds = I, 

Jo 

which implies that rj + (t, x' , I) £ {x 3 = 1} for all t > 0, i.e., that the part of the upper fluid in contact 
with the fixed boundary {x 3 = 1} never flows down from the boundary. It may, however, slip along 
the fixed boundary, since we do not require v + (t, x f , I) ■ ej = for i = 1, 2. A similar result holds for 
rf- at the lower fixed boundary {x 3 = — m}. 

In the steady-state case, the flow map is the identity mapping, rj = Id, so that v = u and q = p. 
This means that the steady-state solution, po, constructed above in Eulerian coordinates is also a 
steady state in Lagrangian coordinates. Now we want to linearize the equations around the steady- 
state solution v = 0, 77 = Id, q = po, for which S- = Id^^o} and A = I. The resulting linearized 
equations read as 

{d t r] = v, 
dtq + podivv = 0, (1.8) 
p d t v + V(p'(p )q) = -gqe 3 - gpoVr] 3 + 2p ajv 2 ei - 2p ujvie 2 . 
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The corresponding linearized jump conditions are 

lv • e 3 ] = and W{po)q} = 0, 

while the boundary conditions are 

v-(t, x , —m) ■ e 3 = v+(t, x' , I) • e 3 = 0. 

1.4 Main results 

Before stating the first result concerning linear problem, we define some terms that will be used 
throughout the paper. For a function / G L 2 (12), we define the horizontal Fourier transform via 

7(6, 6,^3) = I f(x 1 ,x 2 ,x 3 )e-^ +x ^dx 1 dx 2 . 
By the Fubini and Parseval theorems, we have that 

I \f( x )\ 2 dx = ^ f \f(tx 3 )\ 2 dCdx 3 . 

For a function / defined on Q we write /+ for the restriction to + = R 2 x (0,1) and /_ for the 
restriction to il_ = M 2 x (— m, 0). For s£i, define the piecewise Sobolev space of order s by 



H s (n) = {f\u e H s (n+)j„ g H s (n^)} 



(1.9) 



endowed with the norm H/H^s 



given by 



2 



H-{0+) 



+ 



Jj s /Q y For G N we can take the norms to be 



k 

E 

3=0 
k 

E 

3=0 



■xl± 



a + iei 



2^fc-j 



L 2 (I±) 



d£dx 3 



for I_ = (—m,0) and i+ = (0, I). The main difference between the piecewise Sobolev space H s (£l) 
and the ususal Sobolev space is that we do not require functions in the piecewise Sobolev space to 
have weak derivatives across the set {x 3 = 0}. Now, we may state our result on the linear problem 
as follows: 



Theorem 1.1 The linear problem hl.8\) with the corresponding jump and boundary conditions is ill- 
posed in the sense of Hadamard in H k (Q) for every k. More precisly, for any k, j G N with j > k and 
for any Tq > and a > 0, there exists a sequence of solutions {(r) n ,v n ,q n )}™ =1 to hi. 8(1 , satisfying 
the corresponding jump and boundary conditions, so that 



\Vr. 



Hi + \\Vn(P)\\ H i+\\q n (P)\\ Hj < 



n 



(1.10) 



but 



\v n (t)\\ H k > \\r]n(t)\\ H k > a for all t > Tq. 



(1.11) 
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Theorem 1.1 shows discontinuous dependence on the initial data. More precisely, there is a 
sequence of solutions with initial data tending to in H k (Q), but which grow to be arbitrarily large 
in H k (Q). Here we describe the framework of the proof, which is inspired by [3]. First the resulting 
linearized equations have coefficient functions that depend only on the vertical variable, X3 6 (—m, I)- 
This allows us to seek "normal mode" solutions by taking the horizontal Fourier transform of the 
equations and assuming that the solution grows exponentially in time by the factor e^'^*, where 
£ £ M 2 is the horizontal spatial frequency and A(|£|) > 0. We show in Theorem 2.2 that A(|£|) —¥ 00 
in some unbounded domain, the normal modes with a higher spatial frequency grow faster in time, 
thus providing a mechanism for RTI. Indeed, we can form a Fourier synthesis of the normal mode 
solutions constructed for each spatial frequency £ to construct solutions of the linearized equations 
that grow arbitrarily quickly in time, when measured in H k (Q) for any k > 0. Comparing with 
[3], here the main difficulty of constructing this growing solutions lies in building the variational 
structure of the linearized equations, because the natural variational structure breaks down in the 
presence of the rotation term. This difficulty will be circumvented in Section 2 by employing an 
approach that was used first by Guo and Tice [4] to overcome a similar difficulty arising from viscous 
compressible fluids, and later adapted by Jiang, Jiang and Wang [7] in a nontrivial way to construct 
growing mode solutions for viscous incompressible fluids with magnetic field. Due to presence of the 
rotating term, we have to impose stronger restrictions on the parameter s and the spatial frequency 
£ in order to obtain the existence result (|2.2ip and the lower boundedness of A 2 in (|2.24p for the 
rotating case. In addition, the auxiliary function <J>(s) (see [H (3.52)]), constructed by Guo and Tice 
to show (|2.2ip . can not be applied to our case, and we have to construct a new auxiliary function 
F(s) in order to get (|2.21|) . At last, in Section 3 we show a connection between the growth rate of 
solutions to the linearized equations and the eigenvalues A(£), which then gives rise to a uniqueness 
result (see Theorem 3.1). In spite of the uniqueness, the linear problem is still ill-posed in the sense 
of Hadamard in H k (Q) for any k, since solutions do not depend continuously on the initial data. 

With the linear ill-posedness established, we can obtain the ill-posedness of the fully nonlinear 
problem in some sense. We rephrase the nonlinear equations (jl.7l) in a perturbation formulation 
around the steady state, that is, v = 0, r\ = rj" 1 = Id, q = po with A = I and S- = 5+ = Id| X3=0 }. 
Let 

77 = Id + fj, rf 1 = Id — C, v = + v , q = p + a, A = I - B, 

where 

00 

B T = J2(-l) n - l (Dr)) n - 

n=l 

In order to deal with the term h(q) = h(po + a) we introduce the Taylor expansion 

h(p + 0) = h(p ) + ti(p )a + <K, 
where the remainder term is defined by 

/•cr(t,x) 

9K(t,x) = / (a(t,x) — z)h" (po(x) + z)dz 
Jo 

/•po{x)+a(t,x) 

= / (po(x) + a(t,x) - z)h"(z)dz. 

J po(x) 
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Then (|1.7p can be written for fj, v, a as 

d t fj = v, 

d t <J + (po + v)(dwv - tr(BDv)) = 0, 

(1.12) 

dtv + (J - B)V(h'( Po )a + gfj 3 + JH) 

= 2wt; 2 (/ - 5)V(Idi + r?i) - 2w«i(I - 5)V(Id 2 + r/ 2 ), 

where Idj denotes the i-th component of Id, i = 1, 2. We require the compatibility between £ and f/ 
given by 

C = r/o(Id-C). (1.13) 
The jump conditions across the interface are 



J (u+(i,a/,0) - u_(t, 5- (t,x'))) ■ n(t,x',0) = 0, 
\p+(Po + o-+(t,x',0)) =p-(t, P Q +a-(S-(t,x'))), 



(1.14) 

x ,U)) =p-[t,p + a-{i>-{t,x- ))), 

where the slip map (|1.6p is rewritten as 

S- = (H R a - C-) (Mrs + f/+) = Id R 2 + f/+ - £- o (Id R a + ??+). (1.15) 

Finally, we require the boundary condition 

v-(t, x , —m) ■ e3 = v+(t, x' , I) ■ e$ = 0. (1-16) 

We collectively refer to ()1.12|) — (|1.16p as the perturbed problem. To shorten notation, for k > we 
define 



,v,a)(t)\\ Hk = \\v(t)\\m + Ht)\\ Hk + \\<j(t) 



In order to prove the ill-posedness for the perturbed problem by contradiction, we state a definition 
introduced in [3J: 

Definition 1.1 We say that the perturbed problem has property EE(k) for some k > 3 if there exist 
6,to,C > and a function F : [0, 5) — > R + satisfying 

II (77o,uo,cro)|| < S, 
there exist (fj,v,a) G L°°((0, to), H 3 (Q)), such that 

(1) (fj,v,a)(0) = (rjo,v Q ,<T ), 

(2) rj(t) = Id + fj(t) is invertible and rj~ l {t) = Id - Q{t) for0<t< t , 

(3) fj, v, a, C solve the perturbed problem on (0, to) x 

(4) we have the estimate 

sup \\(fj,v,a){t)\\ H 3 < F(\\(fjo,Vo,o- )\\ H k). 
o<t<t 
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Here the EE stands for existence and estimates, i.e. local-in-time existence of solutions for small 
initial data, coupled to L°°(0, to; H 3 (Q)) estimates in terms of H k (Q)-norm of the initial data. If we 
were to add the additional condition that such solutions be unique, then this trio could be considered 
a well-posedness theory for the perturbed problem. 

We can now show that property EE{k) cannot hold for any k > 3. The proof utilizes the Lipschitz 
structure of F to show that property EE(k) would give rise to certain estimates of solutions to the 
linearized equations (11, 8ft that cannot hold in general because of Theorem 1.1. 

Theorem 1.2 The perturbed problem does not have property EE{k) for any k > 3. 

Remark 1.1 Theorem 1.1 and 1.2 show that rotating angular velocity u: can not prevent the linear 
and nonlinear RTI in the sense described in Theorem 1.1 and 1.2, respectively. However, in the 
construction of the normal mode solution to the linearized system in section 2.2, the rotation dose 
have a stabilizing effect on the growth rate X for sufficiently large fixed |£|. In fact, in section 2.2, we 
know that there exist a couple 

(<Po,ipo) e-4= h<P,1>) e L 2 (-m,l) x Hl(-m,l) Po ((p 2 + ip 2 )dx 3 = 1 

such that 

1 f l 

> -A 2 = - / {p'(po)po(^o + |£|<A)) 2 - 2 5 |£|pot/wo + 4^WoM 2 } dx 3 

i r l 

= inf - / {p'(p )po(ip' + |£M 2 - Ht\pv4>V + WWo/A 2 } dx 3 . 

Obviously, ipo ^ 0. On the other hand, we denote the growth rate by Ao for the equations of com- 
pressible inviscid fluids without rotation (i.e., u = 0), that is 

i r l 

> -Xl = inf - / {p'(p )p (^ + |ei^) 2 - 2g\Z\po^} dx 3 . 

{<P,il>)eA 2 J_ m 

Thus we have 

-A 2 > -A§ + 2 I uj 2 pQifl/\ 2 dx 3 > -Xl for cj ^ 0, 

J — m 

which implies X < Xq. 

2 Construction of a growing solution to (11.81) 

We wish to construct a solution to the linearized equations (|1.8p that has a growing if fc -norm for 
any k. We will construct such solutions via Fourier synthesis by first constructing a growing mode 
for a fixed spacial frequency. 

2.1 Growing mode and Fourier transform 

To begin, we make an ansatz 

v(t,x) = w(x)e xt , q(t,x) = q(x)e xt , rj(t,x) = fj(x)e xt , 
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for some A > 0. Substituting this ansatz into (jl,8p . eliminating fj and q by using the first two 
equations, we arrive at the time-invariant system for w = (wi,W2,W3): 

\ 2 p w - V(p'(p )po<livw) = gp divwe 3 - gp Vw 3 + 2p ujw 2 ei - 2p Q ujw l e 2 (2.1) 

with the corresponding jump conditions 

lw 3 ] = and {p (p ) p divw] = 0, 

and the boundary conditions 

w 3 (t,x' ,-m) = w 3 (t,x' ,1) = 0, i'eR 2 . 

Since the jump only occurs in the e 3 direction, we are free to take the horizontal Fourier transform, 
which we denote with either : or J 7 , to reduce to a system of ODEs in X3 for each fixed spacial 
frequency. 

We take the horizontal Fourier transform of wi, w 2 , w 3 in (|2.ip and fix a spacial frequency 
f = (£i,£ 2 ) € M 2 . Define the new unknowns 

tp(x 3 ) = iwi (6>6, x 3 ), 0(x 3 ) = iw 2 (^i,^2,x 3 ), ip(x 3 ) = w 3 (d,C 2 ,x 3 ). 

To write down the equations for (p, 9, ip, we denote ' = d/dx 3 to arrive at the following system of 
ODEs 

' \ 2 p tp + 6[p'(po)/Oo(6^ + 6# + -0')] = tigpoip + 2 Po uj9, 

< A 2 p o + 6[p'(po)po(6^ + 60 + ip')] = &9Po*P - 2 Po ujp, (2.2) 

k AW - W(po)po(£i<p + 60 + #)]' = spo(&p + 60), 

along with the jump conditions 

M = 0, 

W(po)po(Hi<p + 60 + ^)1 = 0, 

and boundary conditions 

^(-m) = V(0 = 0. (2.3) 

We can reduce the complexity of the problem by removing the component 9. For that purpose, 
note that if 99, 0, i\) solve the above equations for £1, 6 an d A, then for any rotation operator 
TZ G 50(2), (ip, 9) := TZ(p, 9) solves the same equations for (6)6) := ^-(6> 6) with tp, A unchanged. 
So, by choosing an appropriate rotation, we may assume without loss of generality that 6 = an d 
6 = \C\ > 0- I n this setting, 9 solves 

A 2 pq9 = —2pouiip. 
Putting this identity into (|2.2p . we arrive at 

>?Po<P + \€\\p'(f>o)Po(\Z\<P + *!>')] = \Z\gp i;-4uj 2 \- 2 p <p, 

(2.4) 

X 2 poip - W(po)po(\Z\<P + 4>')}' = \i\9PW, 
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along with the jump conditions 

(2.5) 

and boundary conditions ()2.3|) . 

In the absence of rotation (oj = 0) and for a fixed spatial frequency £ 7^ 0, the equations (|2,4|) . 
(|2.5p and (|2.3|) can be viewed as an eigenvalue problem with eigenvalue —A 2 . Such a problem has 
a natural variational structure that allows us to construct solutions via the direct method and a 
variational characterization of the eigenvalues via 



where 

1 



1 



and 



1 



^ J — m 

This variational structure was essential to the analysis in [3], where the ill-posedness results for both 
the inviscid linearized problem and the inviscid non-linear problem ( (|1.12|) with u = 0) were shown. 
Unfortunately, when rotation is present, the natural variational structure breaks down. In order to 
circumvent this problem and restore the ability to use the variational method, first we artificially 
remove the dependence of fj2.4|) on A~ 2 by defining s := A~ 2 > 0, and then consider a family (s > 0) 
of the modified problems given by 



-A 2 p </? = l£lb'(Po)Po(l£l</? + V / )] - \€\gpoil> + 4w 2 8po<p, 
-AW = -[p'(Po)po(\t\tp + W - |£|<?/W, 



along with the jump conditions 



M=o, 

W(m)m(K\f + i'')l=o, 



and boundary conditions ()2.3 



2.2 Construction of a solution to ( 12. 6j) 

In the following, we establish the variation framework by formulating constrained minimization. 
Multiplying (|2.6P i and (|2.6p ? by (p and ip, respectively, we add the resulting equations, integrate over 
(— m, I), integrate by parts, and apply the boundary and jump conditions to deduce that 

A 2 f 1 1 f 1 
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Notice that for any fixed £, this is a standard eigenvalue problem for —A 2 . It allows us to use the 
variational method to construct solutions. To this end, we define the energies 

I rl 

E(<P,i>) = 9 / {p'(po)Po('0 / + |£M 2 - 2g\£\p Q ijj(p + 4u! 2 sp <p 2 } dx 3 
^ J—m 

and 

1 f 1 

which are both well-defined on the space L 2 (—m,l) x Hq(— m, I). We define the set 

A={((p,ip) G L 2 (-m,l) x Hl(-m,l) \ J{<f^) = 1}, 

We want to show that the infimum of E(cp, ijj) over the set A can be achieved and is negative, and 
that the minimizer solves the problem (I2.6P , (12. 7p and (|2.3p . Notice that by employing the identity 
— 2ab = (a — b) 2 — (a 2 + b 2 ) and the constraint on J(ip, tp), we may rewrite 



1 



E(<p, 1>) = -g\£\ + - / p'(po)po(rp' + \^) 2 + g\C\po(^P ~ Wdx 3 
+ 2u 2 s ! p^ 2 dx 3 > -g\£\, (<p,i/>) G A. 

J —m 



(2. 



In order to emphasize the dependence on s G (0, oo), we will sometimes write 

E((p,ij}) = E(tp,ip;s) 

and 

p(s) =: -A 2 (s) := M Efafrs). (2.9) 
ipeA 

By inequality (12. 8j) we have 

A < y^. (2.10) 



Proposition 2.1 There exist three constants R\ > 2, Co and C\ depending on the quantities p , 
P±> 9; h m > w ; 80 that 

p(s) < -ColCl + *Ci for any \£\ > R x . (2.11) 

In particular, we have 

p(s) < for any |£| > R\ and s < Cq\(\/C\. 
Proof. Since both E and J are homogeneous of degree 2, it suffices to show that 

mf — 7T- < 0. 
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We assume that <p = — such that the first integrand term in E(tp,ip) vanishes, that is, 

= _ S Mm _ i jy^ dX3 + ^ f^ iX3 ( , 12) 



2 2 J- m p'(m) |{| ! 

where we have used the fact that po solves (jl.4p . Notice that [po] > so that the right-hand side is 
not positive definite. From (|2.12p we see that the rotation term diminishes obviously the growth of 
instability. 

For |£| > 2 we define the test function tp^ E H^—m, I) defined by 

[1 - 30^ , x 3 G [0,/) 
(l + f ) 2 , x 3 g (-m,0). 



Thus, we can estimate that 



and 



for some constants Al, ^2 and A 3 > depending on /9 , pj-, Z, m, g and u, but not on |£|. Making 
use of (|2.14p and (|2.15p . we infer that 

„ , ( < A, - A 5 \£\ + C 1S (2.16) 

for some positive constants C\, A± and ^5 depending on the same parameters, but not on |£|. 

Obviously, there exist a sufficiently large positive constant R\ > and a positive constant Co 
depending on A4 and j4 5 , such that 

A 4 - A 5 |£|/2 < for any |f | > R x . (2.17) 

Inserting (|2.17p into (|2.16p . we obtain 

which implies the desired conclusion. □ 
The key point in the proof argument of Proposition 12. II is to construct a pair (cp,t/j), such that 

ffW 2 (Q) ^ 
2 

which in particular requires that V'(O) 7^ 0. We can show that this property is satisfied by any 
(<p,ip) G A with E((p,ip) < 0. 
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Lemma 2.1 Suppose that (<p,ip) G A satisfies E(<p,ijj) < 0. Then if)(0) ^ 0. 
Proof. A completion of square allows us to write 



2 

9y/P0$ \ . , ,, 9 Po 



V vV(Po) / P^Po) 

Thus, similarly to (j2. 12|) . we can rewrite the energy as 



2 J-m \ S Vp'(Po) 

from which we deduce that if E(ip, if)) < 0, then ip(0) ^ 0. □ 
From now on, we denote 

M := {(£,*) |£ G M 2 , |e| > Ri, s < Coiei/Cl}, 

where the constants R\, Co and C± are from Proposition 12.11 Next, we can show that a minimizer 
exists and that the minimizer satisfies (|2.6|) . (12. 7h and (12.3|) for each (£, s) G M in the same manner 
as in [3]. 

Proposition 2.2 For any fixed (£, s) G At, F achieves its infinimum on A. 

Proof. (|2.8p shows that E is bounded from below on A. Let (</?„, ^/> n ) G „4 be a minimizing sequence. 
Then ip n is bounded in L 2 (—m,l) and 7/? n is bounded in Hq(— m, Z), so up to the extraction of a 
subsequence (p n <p weakly in L 2 , ip n — ^ ^ weakly in fFj, and ip n ^ ip strongly in F 2 . In view of 
the weak lower semi-continuity and the strong F 2 -convergence ip n — > if), we find that 

E((p,ip) < lim inf E((p n ,ip n ) =infF. 

n— >oo A 

All that remains is to show (<p,ip) G A. 

Again by the lower semi-continuity, we see that J(cp,ip) < 1. Suppose by contradiction that 
J((p,if)) < 1. By the homogeneity of J we may find a > 1 so that J(aip,aip) = 1, i.e., we may scale 
up (<p,ip) so that {oap,aip) G A. By Proposition 2.1 we know that inf E < 0, from which we deduce 

that 

E(aip, aif)) = a 2 E(cp, if)) < a 2 inf E < inf F, 

A A 

which is a contradiction since (oup, aif)) G A. Hence J(<p, if)) = 1 so that (<p, if)) £ A. □ 

Proposition 2.3 Let (£, s) G At, and (<p,ip) £ A be the minimizer of E constructed in Proposition 
2.2. Let fj, : = —A 2 := E((p,ip). Then (p,ip satisfy 

\ WW = \i\\p'{Po)Po{\ZW + VO] - \£\9Poi> + 4u 2 sp ip, 

< (2.18) 

[ pp ijj = -\p (p )po(l£lv? + VO]' - If Is/w 

along with the jump conditions \2. 5\) and boundary conditions \2. 3\) . Moreover, the solutions are 
smooth when restricted to either (— m, 0) or (0,/). 
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Proof. Fix ((po, ipo) G L 2 (—m,l) x Hq(— m, I). Define 

j(t, r) = J{ip + tipo + r<p,ip + tip + rip) 
and note that j(0,0) = 1. Moreover, j is smooth, and 

^(0,0)=/ pQ[ipo<p + ihi/)]dx 3 , ^(0,0)=/ p (^ 2 +ip 2 ) = 2. 

•J 711 o — 777. 

Thus, by the inverse function theorem, we can find a C 1 -function r = a(t) in a neighborhood of 0, 
such that <t(0) = and j(t,a(t)) = 1. We may differentiate the last equation to infer that 

^(o,o) + ^(o,oK(o) = o, 



whence, 



<A0) = -^(0,0) = -i j p [^ + ^]dx 3 . 



Since (ip, ip) is a minimizer of E over A, one has 





dt 



E{(p + tipo + cr(t)<p, ip + + cr(t)ip), 
t=o 



which implies that 



J —m 

-S^IPoOOo + o-'(0)tp) + ip(ip + a'{Q)ip)) + 4w 2 sp ¥>(<A) + a' (0)ip)dx 3 . 
The above equation, by rearranging and plugging in the value of ct'(0), can be rewrite as 

p'(po)po(4>' + \£\p)(<Po + l£l<A)) - 5l£|po(Vvo + <P^o) + 4w 2 sp ^o^3 
= P Po[poP + V ; oV']^3, 

J—m 

where the lagrange multiplier (eigenvalue) is p = E((p,ip). 

By making variations with ipo,ipo compactly supported in either (— m, 0) or (0, /), we find that ip 
and ip satisfy the equations (|2.18p in the weak sense in (— m, 0) and (0,/). Standard bootstrapping 
arguments then show that (tp,ip) are in H k (—m,0) (resp. H k (0,l)) for all k > when restricted to 
(— m, 0) (resp. (0,1)), and hence the functions are smooth when restricted to either interval. This 
implies that the equations are also classically satisfied on (— m, 0) and (0, 1). Since (<p,ip) £ H 2 , the 
traces of the functions and their derivatives are well-defined at the endpoints x 3 = —m, and I, and 
it remains to show that the jump conditions are satisfied at 23 = and the boundary conditions 
satisfied at x 3 = —m and I. Making variations with respect to arbitrary ifQ, ip$ £ C£°(—m, I), we find 
that the jump condition 

b>o)po(|^ + ^)J=0. 

must be satisfied. Note that the conditions lip} = and ip(—m) = ip(l) = are satisfied trivially 
since ip G H&(-m,l) ^ C^ 1/2 {-m,l). □ 

The next result establishes continuity property of the eigenvalue p(s) which will be used in finding 
a s with s = 1/A 2 (s) in Theorem 2.1. 
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Proposition 2.4 Let fi : (0,oo) — > M be given by i2.9\) , then fi £ C lo ' c (0, oo), and in particular, 
^eC°(0,oo). 

Proof. Fix a compact interval / = [a, b] C (0, oo), and fix any pair ((fo,ipo) £ A. We may decompose 
E according to 

E(<p,i(>]s)=Eo(<p,il>) + sE 1 (<p t il>), (2.19) 

where 

1 f l 

^ J—m 

Bk(tp,ij}) := 2w 2 / po4> 2 dx 3 > 0. 

The non-negativity of Si implies that E is non-decreasing in s with fixed (93, tp) £ .A. 
Now, by Proposition 2.2, for each s £ (0, 00) we can find a pair (ip s ,ifj s ) G A so that 

E{<p s ,ip s ;s) = inf E(cp,ip; s) = fi(s). 

We deduce from the non- negativity of E\, the minimality of (ip s ,ip s ) and the equality (j2.8[) that 

E(tp ,ipa;b) > E(<f ,ipo;s) > E(ip s ,?p s ;s) > sE^ips,^) - g\£\ 

for all s £ Q, whichi implies that there exists a constant < K = K(a,b,ipo,ipo,g, |£|) < 00, such 
that 

supEi(<p s ,iJj s ) <K. (2.20) 

Let Si £ Q for i = 1,2. Using the minimality of ((p Sl ,ip Sl ) compared to ((p S2 ,ip S2 ), we see that 

H(s{) = E(tp Sl ,ip Sl ;si) < E(ip S2 ,ip S2 ;si). 
Recalling the decomposition (|2.19p . we can bound 

E((p S2 ,ip S2 ;si) < E(cp S2 ,ip S2 ;s 2 ) + |si - s 2 \E 1 (ip S2 ,ip S2 ) 
= /i(s 2 ) + |si - s 2 |£'i((/j S2 ,V' S2 ). 
Putting these two inequalities together and employing (12.20p . we conclude that 

K s i) < M( s 2) + - sal- 
Reversing the role of the indices 1 and 2 in the derivation of the above inequality gives the same 
bound with the indices switched. Therefore, we deduce that 

\fi(si) - n(s 2 )\ < K\si - s 2 |, 
which completes the proof of Proposition 2.4. □ 
To emphasize the dependence on the parameters, we write 

^ = ^(^1^3), ip = ips(\€\,x 3 ), and A = A(|£|, s). 

In view of Propositions 12.3 1 and 12.1] we can state the following existence result of solutions to (|2.6p . 
(H2D and ([23D for each (f , s) £ M. 
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Proposition 2.5 For each (£, s) G A4, i/iere exists a solution (p s (\£\,xz), ip s {\i\, X3) urato A = 
> to i/ie problem \2. 6\) along with the corresponding jump and boundary conditions. More- 
over, these solutions ift s (,\£\)ty an d the solutions are smooth when restricted to either (— m, 0) or 

(0,/). 

In order to prove A — > +00 as |£| — > +00, and give the existence of solutions to the original 
problem (|2.5p . (|2.4p and (|2.3p . we shall further restrict f to satisfy 



|f I > #2 := 2C ( 7 1 7C^ > R lt and C 2 = max{4d, (tfiCo/2) 2 }, 
where Co, C\ and R\ are the constants from Proposition 2.1. Thus we have the following conclusion. 
Theorem 2.1 For each f with |f| > R2 > i/iere exists a s G (0, si) := (0, C0R2/C2), such that 

s = l/A 2 (|f|,s). (2.21) 
Remark 2.1 is easy to check that (f, s) G .M /or any |f| > i?2 o^rf s G (0, si). 
Proof. Recalling — A 2 (s) = /i(s), we define 

F(s) = sA 2 (s) - 1 = -s/i(s) - 1. 
According to Proposition 2.4, -F(s) is continuous. Moreover, we have 

F(0) = -1< 0. (2.22) 

Now if there exists a s G (0, si) such that 

F(s) > 0, (2.23) 

then combining ([2722]) with ([233]) . we can find s G (0,s) such that F(s) = 0, i.e., (pOT]) holds. In 
the following, we verify (|2.23|) . 

According to (|2.1ip and the fact C2 > C%, we see that 

-/x(s) > C R 2 - sd > C R 2 - sC 2 , 

which yields 

F(s) > -C 2 s 2 + C R 2 s - 1 := f(s). 

Recalling the definition of si and R2, it is easy to verify that there exists at least one positive root 
< s < si, such that f(s) = 0. This implies (12331) . □ 

2.3 Construction of a solution to (12.21) 

We may now use Theorem 2.1 to think of s = s(|f |), since we can find s G S so that (|2.2ip holds. As 
such we may also write A = A(|f |) from now on. 

Once Theorem 2.1 is established, we can combine it with Proposition 2.5 to obtain immediately 
a solution to (|2.4|) . and in turn a solution to (|2.2|) for each spacial frequency f with |f | > R2. 
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Theorem 2.2 For £ € K 2 with |£| > i?2, i/iene exists a solution ip = X3), = 0(f, X3), V = 
^(^^3) an d A = A(|f|) > ^° (KM)> such that tp{0) 7^ 0. TTie solution is smooth when restricted to 
(— m, 0) or (0,/), and it is equivariant in f in the sense that if 71 e 50(2) is a rotation operator, 
then 







(Tin K l2 0^ 










^21 ^22 




0(£,X 3 ) 






v° 1 y 




\^x 3 )) 



Moreover, 

A 2 > C3KI - 1 := Cxiei/2 - 1 /or |f | > max {^^, i? 2 l , (2.24) 
where Co, C\ and R2 are the constants from Proposition 2.1 and Theorem 2.1. 

Proof. In view of Theorem 2.1 and Proposition 2.5, there exists a solution (</?(|f|), VKI^Dj $ = 0) to 
(|2.2p for the fixed frequency given by (|£|,0) with |f| > R2. Thus we may find a rotation operator 
U G 50(2) so that K{0 = (|£|,0). Define (pfc s 3 ), 0(f, x 3 )) = ^^(I^Us), 0) and ^,x 3 ) = 
■0(1^1, X3). This gives a solution to (|2.2p for any frequency f with |f| > R2. The equivalence in £ 
follows from the definition. 

We proceed to estimate (|2.1ip . According to (|2.14p . for |£| > i?2 > Ri, one has 

- A 2 < -Co|e| +sd. (2.25) 

The fact that s = 1/A 2 from Theorem 2.1 combined with (|2.25p results in 

A 4 -Co|e|A 2 + d >0, 

which leads to 

A 2 > \ (Coiel + ^C 2 |£| 2 -4C 1 ) > ^C 2 |£| 2 -4C 1 , (2.26) 

or 




(2.27) 



On the other hand, in view of Theorem 2.1, we find that 

A 2 > C0R2/C2 = Cj/2 > y/cl, 

which implies that (|2.27p does not hold. Hence, A 2 has to satisfy (|2,26p . 
For |f I > (Ci + 1)/C , we have 

^C 2 |e| 2 -4C 1 > ^C 2 |cf| 2 -4C 1 |e|+4 = V(Ciiei - 2)2 = d|C| - 2. (2.28) 

By virtue of ([236]) and (f23B|) . the estimate (pT24"|) is obtained by denoting C 3 = C x /2. □ 

Next, we derive an estimate for the i? fc -norm of the solutions (tp,ip) with |f | varying, which will 
be used in the proof of Theorem 2.3 when integrating solutions in a Fourier synthesis. 



Rayleigh-Taylor instability for compressible rotating flows 



19 



Lemma 2.2 Let y?(|£|), V'(ICI) ^ e the solutions to \2. 6\) constructed in Theorem 2.2. Let R<±, Co, C\, 
C3 > be the constants from Theorem 2.2. Fix Ro = max{(Ci + 1)/Cq, 2/C3, R2, 1} G (0,oo). TTien 
/or any £ mf/i |£| > i?o and /or each k > 0, there exists a constant > depending on po, p, g, u, 
I and m, such that 



l// fc (-m,0) + ll^(|£|) llii*(-m,0) + 

Also, there exists a constant Bq > depending on the same parameters, such that for any |£| > ; 



H* ( 0,I) + ll^l)ll^(0,O^^El^- 

j=0 



L 2 (-m,l) 



> Bn. 



(2.30) 



Proof. We begin with the proof of ()2.29|) . For simplicity we will derive an estimate of the i7 fc -norm 
on the interval (0, 1) only. A bound on (— m, 0) follows similarly, and the desired result follows from 
adding the two estimates together. First note that the choice of Rq, when combined with Theorem 
2.2 and (HTTPl) . implies that 



g\t\ < m < 1 - c 3 \z\ = 1 - ^iei - ^ici < - y^i- 



(2.31) 



Also keep in mind that po is smooth on each interval (0, 1) and (—m, 0) and bounded from above 
and below. Throughout the proof we denote by C a generic positive constant depending on the 
appropriate parameters. We proceed by induction on k. For k = the fact that (y(|£|), V'dCl)) £ A. 
implies that there is a constant Aq > depending on the various parameters, so that 

lb(|ei)l|^(o,0 + HV'(|ei)l|^(o,0<^o. 
Suppose now that the bound holds for some k — 1 > 0, i.e., 

k-l 

Mm\H^ W ) + \\m\)\u-H ,D < m-i E 

Define w(\£\) := rf+ipo) Po{ip' + ICl^dCD); where ' = <9 X3 . Recalling /x = -A 2 and s = A -2 , the 
system (|2.6|) implies that 



1 W '(\C\) = - g\£\pov- 

These equations, together with (|2.31|) and the fact that |£| > 1, yield that 
lk(l^l)lliY fc - 1 (o 1 0<^(lb(l^l)ll^- 1 (o ) z) + ll^(lCI)l|^- 1 (o,0, 

K(l^l)[|^- 1( o )0 <c|^| 



(2.32) 



H*-i(o,o + IIVKICDIIfffc-ico.o) ' 



so that ||^(|4|)||^ fe( o,o < C^J =0 jeK- 

On the other hand, the definition of io(|£|) implies that 



H k ~ 1 (0,l) 



< 



w(\a\) 



P'(Po)po 



H k - 1 (0,l) 
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such that ||V ; (|^l)ll_H"fc(o,i) — CEj=o l£P- Returning to the first equation in (|2.32p . we see that 



from which we get 



H k (0,l) 



\i + 4p uj 2 /p 

lei 



fi 1 + 4/9 W 2 / jl 2 



< 



I Ml 



g\\po^{\i\)\\mm + IKI£l)llff fc (o,o 



Thus, making use of (|2.3ip we deduce that 

\u k (Q,l) + 



<C£*=ol£P- ^nce, 



H k (0,l) 



< 



3=0 



for some constant > depending on the parameters, i.e., the bound holds for k. By induction, 
the bound holds for all k > 0. Finally, (pOOj) follows from the fact that {<p(\£\),ip(\£\)) G -4. and p is 
bounded from above and below. □ 

2.4 Fourier synthesis 

Inspired by [3], we will use the Fourier synthesis to build growing solutions to (|1.8|) out of the 
solutions constructed in the previous section for fixed spacial frequency (£l 2 . The solutions will 
be constructed to grow in the piecewise sobolev space of order k, H k , defined by (jl.9p . 

Theorem 2.3 Let Rq < R3 < R4 < 00, where Rq > 1. Let f G C^°(M. 2 ) be a real-valued function, 
such that /(£) = and supp(/) C B(0,R A )/B(0,R 3 ). For £ G R 2 define 

w(£, x 3 ) = -i<p(£, x 3 )ei - i0(£, x 3 )e 2 + ^(£, ^3)e3, 
where ip, 9, ip 7r are solutions to \2.2i) . Writing x' • £ = xi£i + X2C2, we define 



1 



v(t, x) 
q(t,x) 



Air 2 
1 

i^ 2 



(2.33) 



lo O3) 
4vr 2 



/(0(£i¥>(£, ^3) + £ 2 #(£, x 3 ) + d* 3 V>(£, x 3 ))e 



where we have defined A(|£|) = X3) = #(£,£3) = ip(^,x 3 ) = if |£| < -Ro- TTien 77, v, g are 
real-valued solutions to the linearized equations $1.8\) along with the corresponding jump and boundary 
conditions. For every KN, we have the estimate 



IMo)l|jf* + IKo)||^ + Ik 



if* 



<c k { I (l + |£| 2 ) fc+1 |/(£)W < x 



1/2 



(2.34) 



for some constant C k > depending on the parameters po, p, I m, oj and g. Furthermore, for every 
t > we have rj(t),v(t),q(t) G H k (Q) and 



3 tv C4R3 — 1 



lh(0)||^< 11^)11^ < 6*^1117(0)11^, 
jt^nnS^T^^) || Hfc < ||«(t)|| H * < e i v / ^||t;(0)||^, 



(2.35) 



^ etv ^i=T|| g( o)||^ < 11^)1^ <e*^[|g(0)[|^. 
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Proof. For each fixed £ E M 2 , 

V (x,t) = f(Ow^^)e Xmt e lx '<, 
v(x,t) = \mf(Ow(^x 3 )e x ^ t e ix '<, 

q(x, t) = -po(x 3 )f(mMZ, x 3 ) + £ 2 9(Z, x 3 ) + d ? Mi, x^e^e™'* 
gives a solution to (|1.8|) . Since supp(/) C B(0, R±)/B(Q, R 3 ), Lemma 2.2 results in 

sup 

\\dxsw(t, < oo for all k G N. 

£esupp(/) 

These bounds imply that the Fourier synthesis of the solutions given by f|2.33j) is also a solution 
to (jl.8p . The bound ()2.34p follows from Lemma 2.2 with arbitrary k > and the fact that / is 
compactly supported. From (I2.3ip . the estimates (I2.35j) follow. □ 



3 Ill-posedness for the linear problem 

We assume that 77, v, q are the real-valued solutions to (jl.8j) along with the corresponding jump 
and boundary conditions established in Theorem 1.1. Furthermore, suppose that the solutions are 
band-limited at radius R > 0, i.e. that 

(J supp(|??(-,X3)| +\v(-,x 3 )\ + \q(-,x 3 )\) C B(Q,R), 

X3£(—m,l) 

where v denotes the horizontal Fourier transform defined by (jl.3p . We will derive estimates for 
band-limited solutions in terms of R. 

Differentiating the second equation in (|1 .8j) with respect to times t and eliminating the 77 term 
by using the first equation, we obtain 

pod u v - V(p ; (p )p divv) + gp divve 3 + 2p udtv 2 ei - 2p Q ajd t v 1 e2 = (3.36) 

along with the jump and boundary conditions 

ld t v 3 j = 0, and fp \p ) p divvj = 0, 
d t v 3 (t, x , -m) = d t v 3 (t, x', I) = 0. 

The band limited assumption implies that supp(#(-, x 3 )) C B(0,R) for all X3 G (—in, I). The initial 
datum for dtv(0) is given in terms of the initial data q(0) and 77(0) via the second linear equations, 
i.e., 

pdtv(0) = -V(p'(po)g(0)) - gq(0)e 3 - gp Vi] 3 (0) + 2pouv 2 (0)e 1 - 2p ^i(0)e 2 . 

By the standard energy estimate procedure (see [H Section 3]), we can establish three energy 
estimates in the following which ensure the uniqueness result in Theorem 3.1. 

Lemma 3.1 For solutions to \3. 36\) it holds that 

o f (P0\a 12 . P'(po)po 9 2 \ r, f glpo}, ,2 

dt / ir\dtv\ + divv-——v 3 ) =d t —z—\v3\ ■ 
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Lemma 3.2 Let v £ H 1 ^) be band-limited at radius R > and satisfy the boundary conditions 
v 3 (t, x' , — m) = v%(t, x' , I) = 0. Then 

f gfpoj I ,2 p'(po)pq 

Jn 2 P31 ' 2 
where the value of A(R), depending R, is given by (4-2) in J3J Section 4-1]. 

Proposition 3.1 Let v be a solution to $3.36\) along with the corresponding jump and boundary 
conditions that is also band-limited at radius R > 0. Then 

Ht)\\h(u) + IIM*)lli» ( n) < Ce2A(R)t (lN°)l&<n) + 11^(0)11^(0) + Hdiv^(0)||| 2(n) ) 

for a constant C = C(po, I, m,p, g, A(R)) > ; where the value of A{R), depending R, is given by 
(4.2) in [3, Section 4.1} . 

Similar to [3], once we get Proposition 3.1, through constructing the horizontal spatial frequency 
projection operator, we can obtain the uniqueness result. Here we give the proof for the reader's 
convenience. 

Let $ e C^°(R 2 ) be so that < $ < 1, supp($) C B(0, 1), and $(s) = 1 for x 6 B(0, 1/2). For 
R > let <&r be the function defined by &r(x) = Q(x/R). We define the projection operator Pr via 

P R f = T- l ($ R Tf), /GL 2 (0), 

where T- ="■ denotes the horizontal Fourier transform in x' . It is easy to see that Pr satisfies the 
following. 

(1) PrJ is band-limited at radius R. 

(2) Pr is a bounded linear operator on H k (Q) for all k > 0. 

(3) Pr commutes with partial differentiation and multiplication by functions depending only on 
x 3 . 

(4) P R f = for all R > if and only if / = 0. 

Now we are able to prove the uniqueness results on rj, v and q. 

Theorem 3.1 Assume that and (r]2,V2,q2) « r e two solutions to lil.8\) . Then rji = rj2, 

v\ = v 2 and qi = q 2 ■ 

Proof. It suffices to show that solutions to (jl.8p with initial data remain for t > 0. Suppose that 
rj, v are solutions with vanishing initial data. Fix R > and define tjr = Prtj, vr = Prv, q R _ = PrQ. 
The properties of Pr show that 77^, vr, qR are also solutions to (|1.8p but that they are band-limited 
at radius R. Turning to the second order formulation, we find that vr is a solution to (]3.36p with 
initial data vr(0) = dtVR(0) = 0. We may then apply Proposition 3.1 to deduce that 

11^^)11^(0) = \\dtv R (t)\\ L 2 {n) = for all t > 0, 

which shows that T)R(t), VR^t) and q R (t) all vanish for t > 0. Since R is arbitrary, it must hold that 
rj(t), v(t) and q(t) also vanish for t > 0. □ 



divf 



P'(Po 



-v 3 



< 



A 2 {R) 



Po\v 
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The solutions to the linear problem (jl.8p constructed in Theorem 2.3 are sufficiently pathological 
to give rise to a result showing that the solutions depend discontinuously on the initial data. Then, 
in spite of the previous uniqueness result, we obtain that the linear problem is ill-posed in the sense 
of Hadamard. Next, we prove Theorem 1.1. 

Similar to [3j, once we get the pathological solutions to the linear problem (|1.8p constructed 
in Theorem 2.3, we are able to show that the solutions depend discontinuously on the initial data 
described in Theorem 1.1. Here we give the proof for the reader's convenience. 

Proof of Theorem 1.1. Fix j > k > 0, a > 0, To > and let Cj, Rq, Bq, C4 > be the constants 
from Theorem 2.3, Lemma 2.2 and Theorem 2.2 respectively. For each n G N, let R(n) be sufficiently 
large so that R(n) > Rq, CiR(n) — 1 > 1, and 



exp(2T VC^R( 



n 



(1 + (R(n) + l)2)i-fc+l(l + R(n)Y 



> a 2 n 2 C 2 /Bl 



Choose f n G C^°(R 2 ), so that supp(/„) C B(0,R(n) + 1)/B(0, R(n)), f n is real-valued and radial, 
and 

(l + \t\y +1 \fn(\(i\)\ 2 dt = 7 ^. (3-37) 

We may now apply Theorem 2.3 with f n , R3 = R(n), and R4 = R(n) + 1 to find rj n , v n , q n that 
solve (I1.8P with the corresponding jump and boundary conditions, such that r] n , v n , q n G H J (£l) for 
all t > 0. By and the choice of /„ satisfying ([337] , we find that (fTTU]) holds for all n. 

On the other hand, we have 



li 2 (-m,0^ 



Mr )||^> / (i + ICI 2 ) fe l/n(OI 2 e 2ToA(l5l) lk(e,-)li; 

(i + ici 2 ) J+1 i/n(e)i 2 ii^(^-)iii 2( _ m , ^ 



> exp(2T 0V / C^ 



(1 + (R(n) + l)2)i-fc+i 
a 2 n 2 C 2 



>^r^/ (i + iei 2 ) fc i/n(6We = « 2 . 

Bn Jm 2 



J 

Here the first bound is trivial, while the second one follows from supp(/ n ) C B(0,R(n) + 1) and 
A(|£|) > y CiR(n) — 1, and the third one from the choice of R(n) and the lower bound (|2.30[) . Since 
X(\C\) > ^/C 4 R(n) - 1 > 1 on the support of /„, we deduce that 

\\v n {t)\\ 2 H k > \\Vn(t)f Hk > ||i7„(r )||^ for all t > 0. 

from which ([l.lip follows. □ 



4 Proof of Theorem 1.2 



The proof is similar to [3J under necessary modifications. We argue by contradiction. Suppose 
that the perturbed problem has property EE(k) for some k > 3. Let S, to and C > be the constants 
and function provided by the property EE(k). Fix n G N so that n > C, applying Theorem 1.1 with 
this n, To = to/2, k > 3 and a = 1, we can find fj, v, a solving (jl.8p . such that 

\\(fi,v,a)(0)\\ Hk < -, 
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but 

IW*)llff3 > ||»7(t)||ff3 > 1 for t> to/2. (4.1) 

For e > we define t]q = efj(0), Vq = ev(0), and <Tq = ea(0). 

Then for e < Sn, we have ||(?7q, %, ^o)\\H k < So, according to EE(k) there exist rj € ,v £ ,<T £ G 
L°°(0, ioj ^ 3 (^)) that solve the perturbed problem with (^q,Vq,0q) as initial data and that satisfy 
the inequality 

sup \\(fj*,v*,a e )(t)\\ H3 < F(\m,tf Q ,v e Q )\\H>) < Ce\\(fj,v,a)(0)\\ Hk < e. (4.2) 
o<t<t 

Now, defining the rescaled functions ff = fj e /e, v e = v e /e, a e = <r e /e, and rescaling (14, 2ft . we infer 
that 

sup \\(fj e ,v*,a*)(t)\\ H s<l. (4.3) 
0<t<t o 

On the other hand, (fj e ,v e ,a e ) satisfies 
f d t ff = v e , 

d t a € + p Q divv e + e(a e divv e - p tx(B e Dv 6 )) - e 2 (a e tx(B e Dv e )) = 0, . | ( ; 

(dtv e + V(h'(po)a e + gez ■ ff) + V*R e — e(i? <E V(/i / (po)^' e + 9V e ) + B € VW) 
k = 2u{v € 2 e 1 - v\e 2 + e[(v%Vfjl - wfV^) + -B e (ufe 2 - «2 e l)] + e 2 B e (vl\7 fj e 2 - ^Vr/f)} 
with boundary conditions 

a/, —m) • es = v+(t, x' , I) ■ = 0, 

where 



B £ := {I-(I + eDfji)~ l )/e 



and 



9?(t,x) = - (ea e (t,x) - z)h"(p Q (x) + z)dz 

e Jo 

(<j e (i,x) — z)h"(po(x) + ez)dz. 

According to the bound (|4.3|) and the sequential weak-* compactness, we have that up to the 
extraction of a subsequence (which we still denote using only e), 

(fj € ,v e ,a e ) - (fj°,v°,a°) weakly-* in L°°((0, t ), H 3 (n)). 

By lower semicontinuity we know that 

sup m ,v°,a°)(t)\\ H3 <l. (4.5) 

o<t<t 

Note that by construction (rf , v e , a e )(0) = (fj,v,a)(Q). Our goal is to show that (ff,v°,a°) satisfies 
the linearized equations (jl.8p . 
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First of all, we give some estimates on B e and 9v\ We may further assume that e is sufficiently 
small so that 

sup \\ea e (t)\\ L ™ < - inf p (x 3 ) (4.6) 
0<t<t 1 - m < x z< 1 

and e < \/(2K\), where K\ > is the best constant in the inequality ||-Fi?||^4 < i£i||.F||#4 ||-B||#4 
for 3x3 matrix-valued functions F, B. The former condition implies that po + ea € is bounded from 
above and below by positive quantities, whereas the latter guarantees that B e is well-defined and 
uniformly bounded in L°°(0, to; H 2 (Q)) since 



\B e 



n=l 



oo 

H 2 71=1 



oo oo 1 oo 1 

< Ylw-'wmb < E v^mr* < E ^ = 2 - 

n=l n=l n=l 

A straightforward calculation, recalling (j4.5j) and (|4.6p . shows that 

sup ||^l e (t)||jy3 < eK% for some constant K3 > 0. 

o<t<t 

Finally, since Id + erf is invertible, we may define C e y i a (Id + ef) e )~ l = Id — e^ 6 , which implies 
that C = ff o (Id - e( € ). The slip map Si:l 2 xK+^M 2 x {0} is then given by 

Si = Id R2 + efj e + - eC+ o (Id R2 + eff + ). 

The bounds on ff and the equation satisfied by ( e thus imply that 

sup \\St(t) -Id R3 ||i,ao < 2e sup \\fj e (t)\\ Lo ° <2eK 2 sup ||ff (t)||# 3 < 2e# 2 , 
0<<<t 0<t<t 0<t<t 

where K 2 > is the embedding constant for the trace map H 3 (Q) «— >■ L°°(lR 2 x {0}). This bound 
allows us to define the normalized slip map Si := (Si — Id R 2)/e as a well-defined and uniformly 
bounded function in L°°(0, t ; L°°(R 2 x {0})). 

We can now parlay the bounds on fj e , v € , a e , B € , into corresponding bounds on dtfj e , dtv e , dt(J e , 
and some convergence results. Making use of (|4.4jl and the above bounds of the rescaled functions, 
we infer that 

sup \\d t fj e (t)\\ H s = sup \\v e (t)\\ H z < 1, (4.7) 

0<t<t 0<t<t 

lim sup \\d t a e (t) + p divv e (t)\\ H 2 =0, (4.8) 
e ^°o<t<t 

sup \\d t a e (t)\\ H 2 < K 4 , (4.9) 
0<t<*o 

lim sup \\d t v e (t) + V(h'(p ))a £ (t) + #e 3 • ff(t) - 2uv\e x + 2uv\e 2 \ U 2 = 0, (4.10) 
e ^°o<t<t 

sup ||^ e (t)|| H2 <# 5 . (4.11) 

0<t<t Q 

We now turn to some convergence results for the jump conditions. We begin with the second 
equation in ()1.14p . which we expand using the mean- value theorem to get 

P+(Po ) + P+( a +Po + (1 - a%)ea e + )ea% 

= P-(Po)+P-(a-Po + (1 - ° (Id R2 + eSi))«ri o (Id R2 + eSi) ' 
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for functions a± : M + xl 2 -> [0, 1]. From the above bounds on a e and S e we get that 

lim sup \\p' + (a%p+ + (1 - a%)ea%) - p' + (j>t)\\v» = 0, (4.13) 
€ ^°o<t<t 

lim sup \\p'_(atp + (1 - aL)ed± o (Id R2 + eSL)) - MPo = ( 4 - 14 ) 



0<t<t 



and 



lim sup 1 1 at o (Id R 2 + eS e _) - a e _\\ L °c < sup ||Va e (t)|| LO o sup \\eS e (t)\\ L <>° = 0. (4.15) 

£ ^°0<t<i 0<t<t o 0<t<t o 

Since p + (/9q~) = P-{pq ), we may eliminate these terms from the equation (14.12P and divide both sides 
by e; then employing (|4,13p - (j4.15p . we deduce that 

hm sup \\p' + (p+)aX(t)-p^( Po )at(t)\\ L ^ =0. (4.16) 

o<t<t 

For the second equation in (ll,14p we first write the normal at the interface as n e = N e /\N e \ with 

N e = (ei + ed xi fj e + ) x (e 2 + ed X2 fj e + ) 

= e 3 + e(ei x ft^r^. + o^i^ x e 2 ) + e 2 {d Xl fj\ x ^77^) =: e 3 + eiVt 

As e — > we have |Af £ | > 0, so we may rewrite the first equation in (I1.14p as 

{v% - v e o (Id Ra + eS £ )) • (e 3 + e^V e ) = 0. 

Clearly sup ||iV e (t)||ioo is bounded uniformly and 

0<t<i 

lim sup ||77 e o (Id K 2 + eS £ ) — v^_\\l°° < sup ||Z?i) e (t) Hi 00 sup ||e£ e (t)||i,°o = 0, 
e ^°o<t<t 0<t<i 0<t<to 

from which we find that 

lim sup ||e 3 -(t^(t)-tiI(t))|Uoo=0. (4.17) 
e ^°0<t<t 

According to ()4.7|) . (|4.9p and (|4.1ip . by Lions- Aub in Lemma, we see that the set {(fj e , a e , v e )} is 
strongly pre-compact in the spaces L°°(0, to; 12"*"). Therefore, 

{f,a e ,v € ) -> (?? ,a ,w ) strongly in L°°(0, t ; Ft). 

This strong convergence, together with the convergence results (|4.8p and (|4.10p and the equation 
<9(f/ e = v e , implies that 

(d t fj e ,d t v e ,d t a e ) -> [d t fj°,d t v°,d t a ) strongly in L°°(0, t ; F^ (O)), 

and that 

<9 t a° + p divu° = 
k d t v° + V(h'( Po )a° + 5 e 3 • f?°) = 2p un%ei - 2p^v\e 2 . 
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We may pass to the limit in the initial conditions (r] e ,v e ,a € )(0) = (r],v,cr){0) to find that 

(fj°,v°,a°)(0) = (fj,v,a)(0) (4.18) 

as well. 

We now derive the jump and boundary conditions for the limiting functions. The index 11/4 is 
sufficiently large to give the L°°(0, to; L°°)-convergence of (fj e , v e , a e ) when restricted to {X3 = 0}, 
{x3 = — m}, and {x 3 = I}, i.e., the interface and the lower and upper boundaries. Combining this 
with ([USD and ([4TT7]) . we deduce that 

p'+(Po>+ = P'-(P P- on {x 3 = 0}, « - v°_ ) • e 3 = on {x 3 = 0}, 

and 

v\ ■ e% = on {x3 = /}, v°_ ■ e 3 = on {2:3 = — m}. 

Henceforth, (f] ,v ,a ) is a solution to (jl.8p along with the corresponding jump and boundary 
conditions on (0, to) x ^ which satisfies the initial condition (I4.18|) . Thus, in view of the uniqueness 
of the linearized equations Theorem 3.1, we have 

(fj°,v°,a°) = (fj,v,a) on [0,t ) x 0. 

Hence we may chain together inequalities (|4.ip and (14. 5p to get 

2< sup \\(v ,v ,a°)(t)\\ H3 < sup \\(v ,v ,a°)(t)\\ wi <l, 

t /2<t<t 0<t<t 

which is a contradiction. Therefore, the perturbed problem does not have property EE{k) for any 
k > 3. □ 
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